Full Counting Statistics as the Geometry of Two Planes 
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Provided the measuring time is short enough, the full counting statistics (FCS) of the charge 
pumped across a barrier as a result of a series of voltage pulses are shown to be equivalent to 
the geometry of two planes. This formulation leads to the FCS without the need for the usual 
non-equilibrium (Kcldysh) transport theory or the direct computation of the determinant of an 
infinite-dimensional matrix. In the particular case of the application of TV Lorentzian pulses, we 
show the computation of the FCS reduces to the diagonalization of an TV x TV matrix. We also use 
the formulation to compute the core-hole response in the X-ray edge problem and the FCS for a 
square wave pulse-train for the case of low transmission. 
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The experimental and technological importance of a 
full quantum mechanical treatment of the response of 
Fermi systems to a time-dependent perturbation, has 
grown as electronic devices have shrunk. In particular, 
the quantum statistics of charge transfer events induced 
by an applied voltage pulse in simple systems such as 
across a tunneling barrier or along a wire need to be 
understood if the limits to their capacity to encode or 
transmit information are to be found P, 0- They are 
also excellent examples of non-equilibrium systems for 
which in some cases complete theoretical treatments are 
known or are feasible. 

Theoretically, interest has concentrated on the full 
counting statistics (FCS) or generating function, xW-, 
for all moments of the charge distribution There 
are results for the case of an infinite train of periodic 
pulses impinging on a tunneling barrier 0, [|| and, for 
the case of a low transparency barrier with a constant 
bias voltage V, the FCS are also known for non-zero 
temperatures [2J. Particular attention has been paid to 
the case of Lorentzian voltage pulses with 'integer area' 
(# V(t)dt is an integer). If the voltage pulses all 
have the same polarity they generate so-called 'Minimal 
excitation states' (MES). These have been shown to min- 
imise the number of excitations in a one-channel Fermi 
gas required to generate a given signal. They have the 
intriguing property that a train of such pulses retain the 
minimal noise property 0, 01 independent of the separa- 
tion or width of the generating voltage pulses provided 
the pulses are all of the same polarity. When such signals 
impinge on a tunneling barrier their nature is reflected 
in the FCS 0]. 

Here we show that the FCS at zero temperature are 
determined solely by the geometry of two planes. For a 
train of TV Lorentzian pulses these two planes coincide 
except in an TV-dimensional subspace. Computing the 
FCS in this case reduces to the diagonalization of TV x TV 
matrix, which we compute explicitly. The known results 
for x(ty follow directly from this geometric formulation 
circumventing the need to work with a formulation us- 
ing Keldysh Green's functions 0] or solving an auxiliary 
Riemann-Hilbcrt problem to invert singular operators [5( . 
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FIG. 1: (color online) The geometry of the FCS. In the space 
of single-particle states, orbitals occupied at zero temperature 
correspond to eigenstates of h with eigenvalue 1 and define a 
mirror plane shown schematically as the x— axis. Unoccupied 
states at zero temperature (with eigenvalue —1) are in the 
complement of plane, one of these is the y-axis and the re- 
mainder are shown as the z-axis. In the insets the dashed line 
divides states into those above and below the unperturbed 
Fermi level but no other ordering by energy is implied. Ap- 
plication of a voltage pulse, biasing one electrode with respect 
to the other, transforms h — > h — UhU^ . Case (a): (TV„ = 1, 
Nb — see text) The corresponding plane has an added di- 
mension and is shown as the shaded plane. Case (b): (N u — 0, 
Nb = 1) The plane rotates by a/2 about an axis perpendicu- 
lar to h. The axis of rotation and the new occupied state are 
found from [h, e -1 ^'*'] (see text). All other states are rotated 
as a function of time but remain eigenstates of both h and h. 



The mirror planes are defined by the operators h = 
2n — 1 and h = UhU\ where n is the ground state den- 
sity matrix and U is the unitary operator which describes 
the effect of applying a voltage pulse between different 
components of the sytem. Directions corresponding to 
occupied single-particle states are in the corresponding 
mirror plane (eigenvalue +1) while directions correspond- 
ing to unoccupied states arc in its complement and re- 
flected (eigenvalue —1). The planes have the dimension 
of the number of occupied single-particle states. The 
intersection of the two planes (and the intersection of 
the two complements) define the states which are occu- 
pied (or unoccupied) in both the initial and transformed 
states and which therefore contribute nothing to the FCS. 
The remaining states are rotated with respect to one an- 
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other by angles which we denote by The directions 
which are rotated can be identified from the commutator 
[h, e -1 "^*'], which defines a spanning set of states for the 
rotated space. 

We consider as a model system a quantum point con- 
tact between two ideal single channel conductors coupled 
by a tunneling barrier localised around x = and treat 
the effect of a voltage pulse, V(t), applied between the 
two conductors (see Figure [J) . We assume that we can 
neglect the contributions to the noise and other moments 
of the distribution of transferred electrons which scale as 
logio£ where to is the observation time and £ is an ultra- 
violet cutoff set by the Fermi energy. These are present 
even in equilibrium and are associated with the fact that 
the charge number in the left or right electrode is not a 
conserved quantity 0]. The FCS of transmitted charge 
across the point contact (which we take to have reflection 
and transmission probabilities R and T = 1 — R) can be 
obtained from the generating function 



X (A)= £ P„e* 



(1) 



where P n is the probability of n particles being transmit- 
ted (from the left electrode to the right one). 
The principle result of [f| is that 



X(A) = (R + Te iK ' x ) N - 



N,, 



I(l + i?Tsin 2 ^(e iA + e- iA -2; 
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where N u is the number of unidirectional events with 
Ki = ±1 determined by the polarity of the voltage pulse 
and Nb the number of bidirectional events. The geome- 
try for the two types of event is illustrated in Figure Q] 
Consider first the case of a simplified system consisting of 
a single state on either side of the barrier. Suppose now 
that, after taking account of the voltage pulse, the state 
on the left side is occupied and empty on the right. This 
is an example of a single MES 0] or unidirectional [|| 
event and is the case illustrated in Fig[TJa): the /i-planc 
(shown as the x— axis) is actually zero-dimensional for 
this simplified one-state system and the h plane is one- 
dimensional. The FCS for this situation is clearly given 
by x(A) = R + Te lX : P Q = R and P 1 = T in p. 

The bidirectional event is illustrated in Fig [ljb). 
Again consider first the case of a simplified system of 
a single occupied state on cither side of the barrier. If 
the occupied states are equivalent {i.e. they have the 
same energy or the same decomposition by energy), then 
transfer across the barrier is blocked and x(A) = 1. If 
the two states are orthogonal to one another (have dif- 
ferent energies or orthogonal decompositions by energy) 
then each state behaves independently. The FCS for this 
case is a product over that for two independent events 
each of which is of the type illustrated in FigQTa). This 



would give x(A) = (R + Te lX )(R + Te~ iX ). In general, 
the projection of the two occupied states onto one an- 
other is equal to cos 2 ^, where ^ is the angle between 
the two planes in Fig [T] Now, the FCS for this sim- 
plified system is just the weighted average of the two 
previous cases: cos 2 § + sin 2 § (R+ Te lX )(R + T e - iX ) = 



1 + RT sin 2 f(e iA + i 



2). 



The general case of a filled Fermi sea is a combina- 
tion of the two types of events illustrated in Figure [TJ 
This follows from the fact that the operator hh is uni- 
tary and therefore has an orthonormal basis which can 
be contracted from its eigenvectors Q . For positive mean 
transfer from left to right electrode N u > Q, the di- 
mensionality of h is N u greater than that of h. There 
are therefore N u directions with eigenvalue of hh equal 
to -1. The remaining eigenvalues come in pairs (e aj 
with j = 1 . . . Nb) and the 2D plane defined by the cor- 
responding basis vectors, e2j-i X e^j, intersects the two 
mirror-planes defined by h and h in two lines inclined 
at an angle ctj/2. The FCS, x(A), is then given by the 
product over the appropriate factors for the N u and JV& 
events to give ©■ 

The case where Nb and N u are both finite corresponds 
to a train of Lorentzian pulses Q and has attracted a 
lot of attention 0, 0] • Here we show that it is possible 
to identify explicitly the sub-space for which the eigen- 
states of hh have eigenvalues different from unity. The 
computation of the FCS reduces to the diagonalization 
of an N x N matrix where N = N u + 2Nb to find the 
angles oti from its eigenvalues (e lai ). The flux through 
the circuit, (f>(t) = —(e/h) J V(t')dt', can be written in 
this case as 



i<j>(t) 
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t-p* 



(3) 



The poles at p m = t m + %r m have residues A m , which are 
determined by the system of equations 



P n - P r. 



,N. 



(4) 



To within overall phase factors, the effect of the voltage 
pulse is equivalent to a multiplication of left electrode 
states by the factor e 1 ^*' . 

The states which define the space in which the eigen- 
values of hh are different from 1 are \ip m (t)} = l/(t—p m ) 
for m = 1, . . . , N. This follows from considering explic- 
itly the effect of U. All single particle states, ip, satis- 
fying {hh - 1)|V>) = /ie^[/i,e-^]|V>) = 0, do not change 
occupancy under the effect of U (and hence make no con- 
tribution to x(A))- For an arbitrary single-particle state 
tp(i), we find using ([3]) and (0| 



(hh - l)|-0) = si g n ( T m') 



2 A m / A* 



fn 



m,m' 



p"„ 



(5) 
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Here f m is the projection of i(j(t) onto \t/j m (t)): 



fn 



2n 



dt 



{t-p m y 



(6) 



The right-hand side of ([5|) vanishes for all ip(t) orthogonal 
to the set of states {tp n (t), n= 1, . . . N}. 

The FCS for the case of a train of N Lorentzian pulses 
can then be found by considering the effect of hh on 
the subspace of the states {ip n (t)}. We find hh\i{j n ) = 
Y^ m M mn\i>m) , where: 



E 



sign(i 



')A m A* 



(p m -P* m .)iPn ~ P* m >) 



(7) 



The angles a.i which determine the FCS in ([2]) are found 
from the eigenvalues (e Mi ) of the N x N matrix M mn . 

Analytic results for the FCS are known in some simple 
cases [4[: i) The {t„} all have the same sign and ii) N = 
2 with sign(riT2) = — 1. The case of N = 4 has also 
been solved analytically 0]. Case i) corresponds to the 
case of MES. All poles of e 1 ^ are in the same half-plane 
and all N states, \tfj n (t)) = l/(f — p m ) are eigenstates of 
hh with eigenvalue -1. This can be seen from Q using 

sign(r m r I ' n ) = 1, (@J) and by considering de d * 1 - . These 
states define the N non-coplanar directions which are in 
one mirror plane and in the complement of the other. 
The FCS for this case of MES corresponds to N u = N 
unidirectional events in @. 

Case ii) corresponds to a single bidirectional event. 
Here the FCS are just found from the eigenvalues of the 
2x2 matrix M mn . Explicit calculation gives the result 
sin 2 a/2 — | ^~^ | 2 Q. In this case, because the states 
and \1jj2} are orthogonal to each other, it is also par- 
ticularly simple to compute the effect of U on the single- 
particle states explicitly. We work with normalized states 
hM = Cili'i) where d = ^J\pi - p*\. We find that the 
effect of U on the state \4>*) = C\jt — p* is to transform 
it into a linear combination of tpi and ip2 '■ 



i<t>{t) 



Cxijpl -P2) 
C 2 (pi - P2) 



1^2} 



Pi - p 2 

Pi - P2 



01 



(8) 



If \ipl) is in the mirror plane corresponding to h then the 
component of the transformed state proportional to 
is in the complement of the mirror plane and its ampli- 
tude is given by sin ^-e^ 1 — a result previously 
obtained using an operator formalism ■ 

In the general case when Nb is not finite the ge- 
ometric approach simplifies calculations when the trans- 
parency, T, (or the scattering phase shifts in other prob- 
lems) are low, it is possible to deduce the results without 
the need for any direct diagonalization. We illustrate 
this by rederiving the standard result for the orthogonal- 
ity catastrophe or Fermi edge singularity (FES) problem 
and then use the method to generalise results for the 



FCS for a system subjected to a sudden pulse of finite 
duration. 

The FES problem arises in the context of the X-ray 
absorption spectrum of a metal. Here, the form of the 
spectrum is determined by the response of the conduc- 
tion electrons to the sudden appearance of a core hole [T3| 
and has been shown to be the consequence of the orthog- 
onality of the ground states of the Fermi gas with and 
without the local potential due to the core hole ■ The 
core hole spectral function is the Fourier transform with 
respect to tf of the overlap between the state, U(tf)\0), 
where |0) is initial ground state at t = and U(tf) the 
unitary time-development operator of the system. U(tf) 
desribes the effect of the new potential which switches on 
at t = 0. (All states are written in the interaction ba- 
sis with Hq the Hamiltonian of the unperturbed system.) 
This is a one-channel version of the problem we have con- 
sidered. The overlap (0|C/|0) = Y\ k cos a^/2, where the 
e ±ia k are j-jjg eigenvalues of hh. When the rotation an- 
gles afe are small, we can expand the cosine to obtain 
(0|[/|0) ~ exp ^Tv(hh — 1)^. In the usual one-channel 

FES problem, U = e^W where <j){t) = 2(j> Q 6{t)9{t f - t) 
and where (f>o is the phase shift of the core hole poten- 
tial computed on the scattering (unperturbed) states of 
the system. If we work in the time-domain, /i(t, t') = 
-^P( tzti) where P denotes Cauchy principal part, we find 



Tr 



(M-l) . ( e *V«'> -1)*-,.) 



= 2 sin2 '/'o log C*/, 



(9) 



(£ is the usual short-time cutoff of order the Fermi en- 
ergy) and we recover the standard result for this problem: 

(o\u\o) = (^/)-^/- 2 . 

The FCS for a train of step pulses in the low trans- 
parency limit is known to be Poissonian Here 
we will consider a periodic signal with pulses of length 
t p and period S applied to the left electrode with re- 
spect to the right electrode: V(t) = 2(f> n (5(t — nS) — 
5(t — nS — t p )). The geometry of this situation is 
very similar to that of the core hole Green's function 
in the FES problem. The general result @ then gives 
logx(A) = E fe !og[l + Ti? sin 2 a fc /2(e 4A + - 2)]. In 
the limit of TR <C 1, we can expand the logarithm and 



,2 a k 



±Tr(Wi- 1). We find 



Tr (hh - 1 



rn.n 



log 



tp + S(m 



log 



S(m — n — 1 
t p + S(m — r 



S(m — n) + ' 



sin 2 (/> , (10) 



where r is short time cut-off used to characterise the 
delta-function 5(t) = , t 'i + -i\ and N u = for this peri- 



odic train of pulses. In (|10[) . the terms with m = n and 
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m = n + 1 arc divergent for r — > 0, while all other terms 
are convergent. With a measurement time NS we ob- 
tain the usual result log x(A) = 2N^- sin 2 O log l -f (e lX + 
e~ iX - 2). 

Finally we comment on the logarithmic contributions 
to the noise which are present even in the absence of an 
applied voltage pulse and limit possible applications and 
the observability of minimal noise states. When calculat- 
ing the FCS, the system is assumed to have been prepared 
in the ground state of the system with zero transmission 
T = 0. The scattering matrix is taken to be different 
from unity only during the observation time to when the 
transmission becomes non-zero (T > 0) and the number 
of particles on a particular side of the barrier is no longer 
a good quantum number. The geometry of this situa- 
tion is similar to the case of the step pulse in voltage ex- 
cept that the unitary transformation acting on the initial 
states mixes states on both sides of the barrier In 
practice, these logarithmic corrections imply that the ob- 
servation time to for the MES and other effects discussed 
here to be observable above the equilibrium noise must 
be short. For a system with cf = lOmcV and setting 
the equilibrium noise at low temperatures (-^ logepto/h) 
equal to the minimal noise obtained using from a MES 
corresponds to working at frequencies v > 300 MHz. 
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